An exactly solvable model of the BCS-BEC crossover 
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We discuss an integrable model of interacting Fermions in one dimension, tiiat allows an exact 
description of tlie crossover from a BCS- to a Bose-like superfluid. This model bridges the Gaudin- 
Yang model of attractive spin 1/2 Fermions to the Lieb-Liniger model of repulsive Bosons. Using a 
geometric resonance in the one-dimensional scattering length, the inverse coupling constant varies 
from — oo to -|-oo while the system evolves from a BCS-like state through a Tonks gas to a weakly 
interacting Bose gas of dimers. We study the ground state energy, the elementary density and spin 
excitations, and the correlation functions. An experimental realization with cold atoms of such a 
one-dimensional BCS-BEC crossover is proposed. 
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Starting with the historic controversy between Bardeen 
and Schafroth about the proper explanation of super- 
conductivity the crossover from a Bardeen-Cooper- 
SchrieflFer (BCS) superfluid with Cooper pairs, whose size 
is much larger than the inter-particle spacing, to a Bose- 
Einstein Condensate (BEC) of molecules composed of 
tightly bound Fermion pairs has been a basic issue in 
many-body physics Very recently, with the observa- 
tion of molecular condensates near a Feshbach resonance 
in cold atomic Fermi gases Q , this problem has regained 
a lot of attention, since for the first time the crossover can 
be studied in detail. While there is general agreement on 
the fact that the evolution between the BCS and Bose- 
limit is continuous, all existing theories of the crossover 
are approximate In particular, there is no reliable 
description of the most interesting regime near the Fesh- 
bach resonance, where the scattering length diverges. It 
is therefore of considerable interest to have an analyti- 
cally soluble model of the BCS-BEC crossover. This is 
what we will provide in this letter for the particular case 
of one dimension (ID). 

Cold gases allow to realize the BCS-BEC crossover by 
driving a mixture of two spin states with an attractive 
interaction through a Feshbach resonance, beyond which 
a bound state appears in the two particle problem in free 
space. On the BCS side of the crossover, pairs only exist 
in the many-body system due to the Pauli blocking of 
states below the Fermi energy. In one and also in two 
dimensions, the situation is quite different, however, be- 
cause any purely attractive interaction produces a bound 
state already at the two particle level. In fact, contrary 
to the 3D case, its existence is both a necessary and suffi- 
cient condition for a BCS instability ^ . As will be shown 
below, an analog of the 3D crossover can be achieved in 
ID by exploiting a confinement induced resonance (CIR) 
in a tight trap where the effective ID scattering length 
exhibits a resonance caused by the mixing with a closed 
channel bound state in the trap 0]. 

We start by considering the Hamiltonian of the 
Gaudin-Yang (GY) model 7] of a spin 1/2 Fermi gas 



interacting via a short range potential gi5{x): 

i—l * z<j 

where N is the total number of Fermions and m their 
mass. The single dimensionless coupling constant is 
7 = mgi/h^n, where n = N/L is the ID density. For 
attractive interactions, the Hamiltonian describes a 
Luther-Emery liquid When 7 ^ 0~, its ground state 
is a BCS-like state with Cooper pairs, whose size is much 
larger than the average inter-particle spacing 0. The 
strong coupling regime with tightly bound molecules is 
now simply reached by increasing the magnitude of 7. 
The resulting Fermion pairs behave like a hard core Bose 
gas, or equivalently like ID no«-interacting Fermions '10|. 
In this manner, obviously, one never reaches a weakly in- 
teracting BEC as one of the limits of the standard BCS- 
BEC crossover in 3D. This is a consequence of the trivial 
fact that the two-body potential giS{x) has a bound state 
only when gi < 0, but none when gi > 0. As a result, in 
the regime gi > 0, the ground state of the Hamiltonian 
is that of repulsive Fermions and thus has nothing in 
common with a weakly interacting Bose gas of molecules. 

In order to obtain a genuine BCS-BEC crossover in 
ID, it is necessary to exploit a geometric resonance in the 
3D two-body scattering problem with a strong transverse 
confinement. For simplicity, we assume the Fermions to 
be confined in a waveguide with radial frequency uj±/2tt 
and oscillator length a± = y^h/muj±. Moreover, we re- 
strict ourselves to a truly ID situation, where only the 
lowest transverse mode is occupied, requiring huj^ to be 
much larger than the Fermi energy ep- As shown by 
Bergeman et al. [TTI |. in such a situation there is always 
precisely one two-body bound state for the longitudinal 
motion with energy e^, whatever the 3D scattering length 
a. Apart from this bound state, all the scattering prop- 
erties are perfectly described by an effective ID delta 
potential giS{x) with strength 6J: 

gi^2hLu^a{l- Aa/aA_y'^ (2) 
As naively expected, an attractive 3D interaction a < 
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implies a negative value of gi, associated with a bound 
state whose energy et = —mgl/A'h? coincides with the 
exact bound state energy ef, in the limit a/a^ 0. Re- 
markably, gi and thus the binding energy e;, remain finite 
at a Feshbach resonance (a — ±00). Entering the positive 
side a > 0, however, the vanishing of the denominator at 
a±/ a — A ~ 1.0326 ^3 leads to a CIR, where gi jumps 
from —00 to +00 just as in a standard 3D Feshbach reso- 
nance, as discussed recently also in [l^- Now, for gi > 0, 
the short range potential gi5{x) no longer has a bound 
state, however it is still there in the 3D problem. Due 
to the condition hjj_\_ ^ ep (equivalently {na±)'^ <C 1), 
the true bound state energy at CIR, ib = —2huj± is 
much lower in energy than any other relevant scales and 
may formally be taken to minus infinity. 

The ID analog of a Feshbach resonance driven BCS- 
BEC crossover in 3D is therefore described by a modified 
Gaudin-Yang model (^), where for positive gi the repul- 
sive short range potential is supplemented by an addi- 
tional bound state with energy ef, = —00. After crossing 
the CIR at I/7 = 0, the unbreakable Fermion pairs are 
described by a Lieb-Liniger (LL) model jl4j of repulsive 
Bosons. On a formal level, the continuous evolution from 
an attractive Fermi to a repulsive Bose gas in one dimen- 
sion is implicit in the Bethe Ansatz equations of the GY 
and LL models. The ground state for both 7 < and 
7 > may be obtained from the solution of the equations: 



7rp(fc) 
N 



14 

fib 
2 



K 



dq 



ip{q) 



K n 72 + [(fc - q)/nf 

^ dk ,,,h^e 

_^ n Zm 



(3) 



where the quasi-momenta distribution p{k) is normal- 
ized as / dkp{k) = n/2 and Eq is the ground state en- 
ergy. For 7 < 0, they reduce to the GY equations of 
an attractive Fermi gas while for 7 > 0, they reduce to 
the LL equations of a gas of dimers This equiva- 

lence shows that in the whole regime 7 > 0, the den- 
sity of Bosons is hb — 'T-/2, their mass is niB = 2m, 
their coupling 55 is identified with the one of the LL 
model and 7 = mg^ /h'^n. Since one is now dealing with 
dimers, the relation between the coupling constant gB 
and the experimentally accessible parameters a and a±_ 
is no longer given by (0), however. It requires a solution 
of the dimer-dimer scattering problem in the presence 
of a transverse confinement. In the limit a <ti a±, one 
can use the free space result for the effective scattering 
length between point like dimers 035 « 0.6a [l^, lead- 
ing to gB ~ 1.2huj± a. Close to resonance {a±/a > A), a 
plausible ansatz is: 



gB — V2hLuj_a (1 — Aa/a±_) 



(4) 



It follows by assuming that an equation of the form Q 
holds also for dimers and the requirement that the res- 
onances match, while a± is reduced by a factor l/-\/2- 
Note that Q corresponds to an effective dimer scatter- 
ing length asB = a/ V2. 
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FIG. 1: Ground state effective energy per particle Eg^ ^ /N = 
Eq/N — eb/2 [in units of the Fermi energy for the non- 
interacting gas ep = TT^h^n'^ /8m] as a function of I/7. 



The qualitative physics of the modified Gaudin-Yang 
model is now simple to understand: when 1/7^ —00 
(BCS limit), the system consists of weakly bound Cooper 
pairs. The associated excitation gap is related to the 
spin sector, i.e. there is a finite gap between the singlet 
ground state and the first triplet excited state. In ad- 
dition, there are gapless density fluctuations describing 
the Bogoliubov- Anderson mode of a neutral superfluid. 
At resonance, when 1/7 = 0, the system is a Tonks gas 
of tightly bound dimers ^3 ■ It still exhibits sound 
modes with a linear spectrum, however the spin sector 
has disappeared because the spin gap diverges. On the 
positive side of the resonance, the system is an interact- 
ing Bose gas of tightly bound molecules. Its excitations 
are the standard Bogoliubov sound modes, whose veloc- 
ity vanishes asymptotically in the weak coupling limit 
1/7 ^ -foo (BEG limit). 

For quantitative results, we start with the ground state 
where the number of up and down spins is identical. Sub- 
tracting the bound state contribution, the relevant finite 
quantity to consider is Eq^-^ = Eq — Neb/2, which we 
call the effective ground state energy. This quantity is 
obtained by numerically solving the Bethe Ansatz equa- 
tions and is plotted as a function of I/7 in Figure 
[T] Its asymptotic behavior in the BCS |3| and BEC-limit 
[Tij l and near the CIR is: 
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im is the Fermi energy for the non- 



where ep 

interacting gas. We note that on resonance E^ 

E'o"h = 0-)/4. 

We now discuss the low energy elementary excitations. 
First consider the density excitations (the "charge sec- 
tor" ) which are gapless phonons with dispersion relation 
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FIG. 2: Sound velocity Vc [in units of the Fermi velocity for 
the non-interacting gas vp = nhn/2m] as a function of I/7. 



FIG. 3: Spin gap A [in units of the Fermi energy for the 
non-interacting gas ep — 7r^fi^n^/8m] as a function of I/7. 



Lu = Vc\k\, when \k\ 0. An effective theory is provided 
by a Luttinger hquid description with charge velocity Vc 
and correlation functions described by the Luttinger pa- 
rameter Kc jl2]. When 7 < 0, Kc'^^ describes fermionic 

correlation functions and is obtained from Kc — vf/vq 
[Tsj j. where vp = TThn/2m is the Fermi velocity of the 

non-interacting gas. When 7 > 0, Kc describes bosonic 
correlation functions and is given by Kc^^ = vf/'2vc [l8l| . 
In order to determine these parameters, it is therefore 
enough to extract the sound velocity from the exact so- 
lution. This is done in the usual way by calculating the 
compressibility from the ground state energy 0|. Its 
asymptotic behavior in the BCS [ij and BEC-limit 
and near the CIR is: 
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Figure 121 shows Vc as a function of I/7. Note that is 
always larger than one, which implies that the system is 
a ID superfluid, as discussed, e.g., in poj j. 

The low energy spin excitations (the "spin sector" ) are 
described by a sine-Gordon model [13 with coupling pa- 
rameter (3. A recent renormalization group analysis of 
the sine-Gordon model 21, 22] shows that when 7 < 
the system is driven to the strong coupling fixed point 
= 47r, i.e., the spin Luttinger parameter Kg |l3| is 
equal to 1/2. There it becomes equivalent to a non- 
interacting massive Thirring model |23j, i.e., to a gas of 
massive relativistic fermions. The latter are called mas- 
sive spinous and may be interpreted as a quantum soliton 
of the sine-Gordon model. They obey the relativistic dis- 
persion relation: 



(7) 



when \k\ — > 0. Therefore, the low energy part of the spin 
sector is fully described by the spin velocity Vs and the 



spin gap A or, equivalently, by the mass of the spinon 
■nis = A/(2-y^). The spin gap is also defined as the en- 
ergy difference between the singlet ground state and the 
first triplet excited state. The spin velocity and the spin 
gap can be extracted from the spinon dispersion relation 
computed from the exact solution of the GY model. The 
spin gap has the following limiting behavior: 
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where the BCS limit was already obtained in When 
7 > 0, the spin gap is infinite in our model. Its be- 
havior as a function of I/7 is plotted in Figure |3| When 
1/7 < —1, it behaves similarly to the gap computed from 
the BCS mean-field theory Abcs ~ e_F exp (— 7r^/2|7|) 
and the dispersion relation is reminiscent of that of 
Bogoliubov quasiparticles, provided Abcs is identified 
with A/2 9]. Near resonance, the spin gap is equal 
to the modulus of the two-body bound state energy 
jefjj = 27^ei?/7r^. The smooth crossover from Cooper 
pairs to molecules occurs when the size of a molecule 
is of the order of the average distance between particles 
n~^, which happens for 7 —2. 

The spin velocity can be computed with the bosoniza- 
tion approach in the weak coupling limit 17] and from 
the Bethe Ansatz equations in the strong coupling limit 
0. It turns out to be given by: 
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and is plotted as a function of 1 /j in Figure ^ Spin 
correlation functions can be obtained from the knowledge 
of the spin gap A and of the charge Luttinger parameters 
Kc and v,,, as shown in 24]. 

The above scenario for a BCS-BEC crossover can be 
realized in an experiment with ultra-cold gases confined 
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FIG. 4: Spin velocity Vs [in units ol the Fermi velocity for the 
non-interacting gas vp ~ nhn/2m] as a function of I/7. 



in a quasi-lD geometry, e.g., ^Li in an array of ID tubes 
created with optical lattices 25] or on an atom chip [2^ . 
In order to reach a ID regime, the quantum of radial os- 
cillation hu!± has to be much larger than (i) the quantum 
of axial oscillation huj^^, (ii) the Fermi energy ep and (iii) 
the thermal energy fc^T. Having a degenerate system 
imposes that the latter is also smaller than the Fermi en- 
ergy. Note that condition (ii) is required for the validity 
of the model, as previously mentioned. Taking, e.g., the 
following values for radial uj±/2tt ~ 100 kHz and axial 
a;||/27r ^ 200 Hz trap frequencies, number of particles 
(per tube) iV ~ 100 |2|| and temperature T ~ 50 nK 
allows to fulfill the previous requirements. To study the 
crossover in such a setting, the dimensionless coupling 
constant 7 needs to be tuned through a Feshbach reso- 
nance 0, giving rise to a CIR as discussed above. The 



experimental characterization of the different regimes can 
be done, for example, by measuring the axial collective 
modes in the trap. The ratio of the frequencies of the 
breathing and dipole modes [2^ is 2 in both the BCS 
and CIR limit 0, and V3 in the EEC hmit In 
addition, extending an available RF spectroscopy tech- 
nique used to measure the molecular binding energy |28j| , 
it should be possible to extract information on the gap. 

In conclusion, we have presented an exactly solvable 
model for a BCS-BEC crossover in ID. The model may 
be realized experimentally by using the combination of 
a standard Feshbach- with a confinement-induced scat- 
tering resonance. In contrast to the 3D case, the exact 
solution allows to make quantitative predictions in the 
whole crossover regime. Our results show explicitly that 
the evolution from the BCS- to the BEC-limit is contin- 
uous, as expected on general grounds. Of course some of 
our results are specific for one dimensional systems, e.g., 
the pairs are unbreakable on the positive side of the reso- 
nance, which simplifies the description there to a purely 
bosonic model. Regarding the behavior of the sound ve- 
locity, however, the situation should be qualitatively sim- 
ilar to the 3D case. As such, our model provides a novel 
and experimentally accessible tool to address one of the 
long standing basic problems in many-body physics. 

We thank W. Rantner and S. Cerrito for useful discus- 
sions and for carefully reading the manuscript. We also 
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Note added. After this work was completed, we became 
aware of a related work of I.V. Tokatly [2^ . 
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